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Research  on  Minimim-Fuel  Fixed  Thrust  Transfers  Between  Coplanar 
Orbits 


Wang  Xiaojun,  Wu  Delong,  Yu  Menglun 

(Beijing  Aeronautics  System  Engineering  Research  Institute, 
Beijing,  100076) 


Abstract:  This  paper  studies  a  minimum-fuel  thrust 
direction  control  strategy  for  a  primary  insertion  from 
a  circular  (or  elliptical)  orbit  to  a  coplanar  (or 
elliptical)  orbit  with  a  fixed  thrust.  Problems  of 
this  kind  can  be  incorporated  into  the  double-bound 
value  problem,  which  is  discussed  with  reference  to  the 
selection  of  initial  free  values  and  the  burn  method  is 
used  for  an  iterative  solution.  Also  in  this  paper, 
optimal  transfer  from  a  parking  orbit  to  a 
geosynchronous  transfer  orbit,  and  between  two  near- 
earth  circular  orbits  are  calculated,  with  satisfactory 
results . 

Key  words:  orbit  transfer,  minimum  fuel  consumption, 
principle  of  the  maximum,  optimal  control 


1 .  Introduction 


Spacecraft  are  generally  required  to  accomplish  one  or 
several  dynamic  orbit  transfers  before  entering  a  target  orbit 
and  therefore  it  is  vitally  important  to  design  an  orbit  transfer 
process  which  spares  propellant  and  reduces  the  spacecraft  mass. 
The  conventional  design  of  orbit  transfer  is  based  on  the  Hohmann 
optimal  orbit  transfer  theory  under  the  fundamental  assumption 
that  velocity  variation  is  completed  in  an  instant,  i.e.  an 
impulsive  orbit  transfer.  This  is  fairly  accurate  with  a  larger 
thrust.  However,  given  the  diversification  of  flight  missions, 
the  requirements  for  overload  becomes  more  and  more  rigorous  and 
the  amount  of  thrust  is  restricted.  Under  such  a  scenario,  more 
and  more  attention  has  been  given  to  limited  or  small  thrust 
multiple  orbit  transfers  during  extended  upper  level  space  flight 
missions  or  interstellar  probes  [1].  When  propelled  with  a  small 
thrust,  a  spacecraft  must  fly  dynamically  long  distances  in  space 
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•and  may  suffer  from  velocity  losses  caused  by  gravity.  In  this 
case,  a  major  concern  is  how  to  control  the  variation  of  the 
thrust  vector  so  as  to  reduce  the  velocity  losses  and  much 
research  has  been  done  on  this  aspect  [2,3]. 

The  variational  method  and  the  principle  of  the  maximum  are 
effective  methods  to  solve,  respectively,  the  problem  of 
liberating  control  of  constraints  and  the  problem  of  optimizing 
the  control  constraints .  An  analytical  solution  can  be  obtained 
for  a  simple  situation,  yet  a  complex  multivariable  problem  such 
as  space  track  optimization  requires  an  iterative  solution  of 
double-bound  values  [3-7],  The  major  difficulty  in  this  solution 
is  how  to  select  the  initial  values  of  those  conjugate  variables 
that  are  not  physically  meaningful,  while  the  equation  is 
extremely  nonlinear  and  its  convergence  depends  on  the  accurate 
selection  of  the  initial  values  [3,5,7].  Thus,  it  is  very 
important  to  transform  the  initial  free  values  into  physically 
meaningful  variables  and  adopt  some  efficient  measures  in 
calculations  based  on  the  possible  optimal  results . 

In  this  paper  the  authors  derive  the  optimal  process  and 
analyze  in  detail  the  bound-value  conditions  for  various  initial 
orbits  and  target  orbits  based  on  the  different  characteristics 
of  circular  orbits  and  elliptical  orbits.  Also  discussed  is  how 
to  select  initial  free-value  variables  using  the  initial 
integration  descending  order.  Finally  the  calculations  are 
carried  out  for  the  fixed  thrust  optimal  transfer  from  a  parking 
orbit  to  a  geosynchronous  transfer  orbit,  and  between  two 
circular  orbits,  together  with  an  analysis  of  the  optimal 
results. 

2.  Theoretical  Analysis 

By  using  a  reference  radius  r^f/  reference  time 
(gref  is  gravitational  acceleration  at  r^ef)  and  initial  mass  mo  to 
make  the  equation  non-dimensional,  a  non-dimensional  dynamic 
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equation  in  the  inverse  ratio  squared  gravitational  field  can  be 
derived  [ 2 ] : 

—  =  Vsin/ 

dr 

d6  _  V^cos/ 

dt  r  ( 1 ) 

dV  _  sin/  ^  cosa 
dr  ~  m 

^  _  j  1_\  cos/  „  sing 

dr  \  r  /  V  mV 

dm  ^1 
dr  ~ 

where  r,  0,  V,  y ,  m,  respectively,  are  the  geocentric  distance, 
polar  angle,  velocity,  flight  track  angle  and  mass,  which  are 
referred  to  as  state  variables;  is  the  engine  gas  injection 
speed  which  is  taken  as  a  constant;  G=T/mogref,  T  is  engine 
thrust;  Earth's  radius  is  usually  selected  as  then  G  is  the 

initial  thrust-weight  ratio;  a  is  the  angle  of  attack. 

The  goal  of  optimization  is  minimum  fuel  consumption  itip, 
which  is  equivalent  to  the  maximum  orbital  insertion  mass  m(Tf); 
the  angle  of  attack  a  is  a  control  parameter  to  be  optimized. 


By  introducing  the  Hamiltonian  function  in  Eq.  (1): 


H  =  /IV'siny 


A,  j  - 


ycos7 


1  i  cos/  ^  sing"]  ,  ^  J_ 
7 !  ^  mV  J  y. 


-t.  -  are  the  conjugate  variables  corresponding  to  the 

state  variables,  which  must  obey  the  following  differential 
equation  along  the  optimal  flight  track  in  accordance  with  the 
optimization  theory  [2]: 


dXr  ,  Vcos/ 


2  1 cos7 


=  —  A,.sin7 


COST  _  ,  [(  JL  ,  _L 


=  —  A^Vcos/  T  A,  -^siny  +  A,r  ^cos7  +  A,|  - i 
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dr  m* 


Xircosa  +  Xr ' 


V 


From  the  optimal  principle: 


an 

da 


is  derived: 


a  =  tan  ’ 


jL 

a,-F 


(3) 


(4) 


The  Hamiltonian  function  H  does  not  contain  distinct  time  x 
and  is  an  autonomous  system  [2];  the  terminal  time  Xf  is  free  and 
the  performance  index  and  terminal  constraint  do  not  contain 
distinct  Xf.  Therefore: 

(=> 

From  Eq.  (3)^-0,  (6) 

=  C, 


Based  on  Eqs.  (1),  (3)  and  (4),  if  the  initial  values  of 
state  variables  and  conjugate  variables  are  known,  the  solution 
to  the  problem  can  be  derived  through  numerical  integration. 

3.  Analysis  of  Bound  Value  Conditions 

The  orbital  constraint  mode  is  given  by  orbit  parameters, 
including  eccentricity  e,  energy  C  and  argument  of  perigee  e. 
Their  relationship  with  state  variables  is  [2]: 

e  =  y  1  -r  v(v  —  2)cos-7 ,  ^  =  \ - w  =  d  —  tan^' 

2  r  p  ~  r 

where  p  _  =  rF- 


With  different  initial  orbits  and  target  orbits,  the 
corresponding  bound  value  problems  are  also  different. 
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3.1  The  Case  When  the  Target  Orbit  is  an  Elliptical  Orbit 


The  parameters  of  target  orbit  are  Qf,  Ctr  Wf  The  state 
variables  should  satisfy  the  constraint  equation  in  the  terminal 
time  as  follows: 


i(r/)  =  V  1  v(tj)(u(T.)  —  2)ccs-7(r/)  —  e,  =  0 


'I'M  -  -  f/  =  0 


=  *9(r/)  —  tan 


/>(r/)tany(rf) 


pit/)  —  r{T,)  ’ 


—  ca,  =  0 


(7) 


then  the  terminal  transverse  conditions  are: 


3^1 


3^2 

driVf) 


3^2 

driZf)  ' 


3^,  ,  3^3 

Av(r/)  -  aV^(r/) 


3^3 


A^(r/)  =  Us 


^  ,  ^3 

WUT)  ^  snvf) 


From  the  above,  the  following  equation  can  be  derived: 


A,(r/) 


1 


r^(r/)V(r/) 

1 


r-(r/)tan7(r/) 


Av(^/)  “h 

/ 

1 


1 


1 


riVfW^iTj) 


Ar(r/^)  =  0 


(8) 


Eqs.  (7)  and  (8)  contain  four  constraining  equations. 

3.1.1  The  Case  When  the  Initial  Orbit  is  an  Elliptical  Orbit 

The  orbit  parameters  are  eo,  Co/  “o*  If  the  orbital  transfer 
location  is  not  restricted,  they  should  meet  the  following 

f 

constraints  in  the  initial  time: 
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^.,(ro)  =  Vl  +  y(ro)(y(ro)  -  2)cos'7(r„)  -  eo  =  0 


^z(To)  =  Y^nr„)  -  -  0 

_,  ^(r.Otan/Cro) 
^,(ro)  =  0(ro)  -  tan- 

■  =  m(To)  —  1  =  0 


OJo  =  0 


(9) 


The  transverse  conditions  at  the  initial  laoinent  arei 


A.(ro)  ^ 

= 

dritf)) 

"^-a^(ro)  ' 

Av(ro) 

av(ro) 

'  ^'aV(ro) 

Ar(rc) 

^(ro) 

'  a/CrJ  ■ 

,  A«(r,,)  =  fii 


a^>. 


■  aF(ro) 

A„(r„)  =  /i, 


Then  the  relationship  similar  to  Eq.  (8)  is  derived  as: 


A>(rr,) 


—  r(rc)tan/(r(,)A,(r„) 


tan/Cr, ) 
r(r,j>V(rcJ 


Ar(r,) 


1 

r(ro)V-(r„) 


iA,(r,,) 


(10) 


Since  (3)  differential  equation  of  conjugate  variables  is 
linearly  homogenous,  the  initial  value  of  one  conjugate  variable 
can  be  determined  randomly  (with  the  required  symbol)  [4],  which 
is  taken  as: 

(11) 

=  1.0 


By  substituting  Eqs.  (10)  and  (11)  in  (5)  and  forming 
simultaneous  a  set  if  equations  with  Eq.  (4),  then: 


(12) 


Ar(rc,)  = 


V(ro) 

V. 


sina(ro)»  AvCr,)  =  ^cosffCrj) 
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The  four  equations  (9-12)  contain  eight  equations  and  eleven 
unknown  variables  (including  a),  of  which  three  variables  cannot 
be  determined,  supposedly  6(to),  a(T:o),  ^.^(To),  referred  to  as 
initial  free  values.  There  are  four  constraining  conditions  at 
termination:  one  of  them  is  the  integration  termination  condition 
(shut-down  condition),  while  the  other  three  constraints 
correspond  to  the  three  initial  free  values,  formulating  a 
double-bound  value  problem. 

3.1.2  The  Case  When  the  Initial  Orbit  is  a  Circular  Orbit 


The  orbit  parameter  is  then  the  constraint  is: 

^i(ro)  =  riTo)  -i-  ~  =  0,  ^2(^0)  =  V(r^)  —  V—  2?o  ==  0 

^  25o 

.^3(^0)  =  ^(^0)  =  0,  =  W2(r„)  —  1  =  0 

and  the  transverse  condition  is : 

^(^0)  =  0 


(13) 


(14) 


Eqs .  (11-14)  contain  a  total  of  eight  equations 
as  above,  and  the  free  initial  values  are  still 
e(ro),  a(To)  and  X^(To) 


3 . 2  The  Case  When  the  Target  Orbit  is  a  Circular  Orbit 

The  orbit  parameter  is  Cf/  and  terminal  constraint  is: 

'f',  <r,>  .  Kr,)  +  -  0.  -V<r,>-  4^,  -  0.  {-.(r,)  _  y^,,^  .  „ 


The  transverse  condition  is:  =  0 

From  Eq.  (6),  =  0,  -  o 

For  the  case  when  the  target  orbit  is  a  circular  orbit. 
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.  there  are  only  three  constraint  conditions  in  Eq.  (15)  at  the 
termination . 


3.2.1  The  Case  When  the  Initial  Orbit  is  an  Elliptical  Orbit 

The  initial  constraint  equation  totally  conforms  to  Eq.  (9) 
By  substituting  Eqs.  (16)  in  (10): 

~  Kr,)tan/'(r„)(  ^  ~  ^  ^ 


Thus,  (9),  (11),  (12),  (16)  and  (17)  contain  nine  equations 
and  two  initial  free  values  selected  as  6(To)  and  a(To). 


3.2.2  The  Case  When  the  Initial  Orbit  is  a  Circular  Orbit 

When  both  the  initial  orbit  and  target  orbit  are  circular 
orbits,  the  location  0  has  no  effect  on  the  evaluation  of  other 
variables  at  the  initial  point  and  termination  and  A0(t)=O. 

Thus,  0  and  Xq  can  be  omitted  with  nine  variables  remaining; 
Eqs.  (11),  (12)  and  (13)  contain  seven  equations  with  two  initial 
free  values  selected  as  a(To)  and  X^{Xq)  . 


4 .  Sample  Calculation  and  Analysis  of  the  Result 

Sample  Calculation  1  Transfer  from  a  parking  orbit  (a  circular 
orbit  200km  in  height)  to  a  synchronous  transfer  orbit  (apogee 
ha=35786km,  perigee  hp=200km,  0=180°).  Ve=4500m/s. 
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Fig.  1  Optimal  Angle  of  Attack  Variation  Law 
Key:  (1)  Angle  of  attack;  (2)  Time 


Fig.  2  Pitch  Angle  variation  Law 
Key:  (1)  Pitch  angle;  (2)  Time 

Figure  1  shows  the  optimal  variation  law  of  attack  angle  a, 
which  is  close  to  a  linear  variation  with  extremely  large 
variation  amplitude  of  a  when  the  thrust  is  small.  Fig.  2 
displays  the  variation  law  of  pitch  angle  (p(t)  in  inertial  space 
when  G  varies  ((p=  0+9O°-a-Y);  «P  varies  in  a  smaller  amplitude 
despite  the  large  variation  of  a  with  a  small  thrust.  Generally 
speaking,  <P(t)  is  designed  to  be  in  a  linear  variation  [8]  in 
©ngineering .  Obviously,  with  a  small  thrust,  the  number  of 
segments  should  be  increased  so  as  to  obtain  an  optimal  control 
law. 


9 


lOOC.  2^  . 

5 


i 

400  *3 

I  G  -  0.2 

200  ^ 

i 

140  160 


0-3.5 


2 

««(*) 


180  200  220 


Fig.  3  Variation  of  Flight  Altitude 
Key:  (1)  Flight  altitude;  (2)  Polar  angle 
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Fig.  4.  Relationship  between  Required 
Velocity  Increment  and  Thrust-Weight  Ratio 

Key:  (1)  Velocity  increment; 

(2)  Initial  thrust-weight  ratio 


Figure  3  shows  the  spacecraft  altitude  variation  law  along 
the  optimal  flight  track:  it  descends  at  first  and  then  ascends 
to  enter  the  target  orbit;  when  G=0.2,  the  minimum  altitude  is 
176km,  and  if  G  decreases  further,  the  minimum  altitude  of  the 
spacecraft  will  decrease  and  even  intersect  with  the  ground.  In 
consideration  of  actual  needs,  such  as  the  need  for  avoiding  a 
dense  atmosphere,  it  is  required  to  restrict  the  minimum  altitude 
(r(t)>r^n)*  This  is  a  problem  of  state  variable  constraints 
along  the  flight  track  (path  constraint),  which  is  expected  to  be 
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.studied  further. 


Figure  4  shows  the  relationship  between  the  velocity 
increment  needed  for  the  spacecraft  to  enter  the  target  orbit  and 
the  initial  thrust-weight  ratio.  The  straight  line  in  the  figure 
corresponds  to  an  ideal  situation  of  impulsive  orbital  insertion, 
while  with  limited  thrust,  Av  is  always  overstated.  This 
difference  is  just  the  velocity  loss  caused  by  gravity — the 
smaller  the  thrust,  the  larger  the  loss  can  be. 

Sample  Calculation  2  Transfer  from  a  circular  orbit  (with 
height  200km)  to  an  adjacent  circular  orbit  (500km  in  height). 
Vg=2800m/s . 

In  launching  some  satellites,  the  launch  vehicle  is  required 
to  send  them  onto  a  higher  orbit.  In  such  cases,  the  mission 
undertaken  by  the  final  stage  of  the  rocket  is  similar  to  a  fixed 
thrust  primary  entry  between  two  circular  orbits.  More  commonly, 
this  kind  of  activity  can  be  referred  to  as  an  optimal  transfer 
between  non-intersecting  orbits.  This  transfer,  while 
accomplished  in  an  impulsive  mode,  requires  an  additional 
transitional  orbit  for  a  secondary  orbit  transfer. 


J60  y .  -  - 

]  «A(°)  1 

G^o  j 

G  =  0  5 

270  j  • 

225  i  :  ;  / 

i  ;  I 

“I  ‘ 

I  80  J  ■  -  r  ; 


1  if,  H 


....  2 

?Oil.i  I'.oii  .itKiO 


Fig.  5.  Optimal  Angle  of  Attack  Variation  Law 
Key:  (1)  Angle  of  attack;  (2)  Time 
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Fig.  6.  Relationship  between  Required 
Velocity  Increment  and  Thrust-Weight  Ratio 

Key:  (1)  Velocity  increment;  (2)  Primary 
orbital  insertion;  (3)  Secondary  orbital  transfer* 
(4)  Initial  thrust-weight  ratio  ' 


Fig.  5  shows  an  optimal  variation  law  of  attack  angle  a,  by 
which  the  spacecraft  is  flying  backward  during  the  period  of  time 
90°<a<270°.  a  varies  dramatically  near  180°  and  it  varies  almost 
abruptly  with  a  larger  thrust.  It  can  be  foreseen  that  if  a  is 
restricted  (I  a(t)l  ,  it  will  vary  abruptly  here.  Fig.  6 

shows  the  relationship  between  the  velocity  increment  needed  for 
an  orbital  transfer  and  G.  During  a  primary  orbital  insertion, 

Av  will  increase  rapidly  with  increase  in  G  and  therefore,  a 
small  thrust  propulsion  is  suitable  for  this  occasion.  For 
comparison,  a  transitional  orbit  (h^  is  500km,  hp  is  200km,  o  is 
180°)  is  utilized  for  studying  the  secondary  optimal  orbital 
transfer,  for  which  a  velocity  increment  must  be  greatly  reduced 
as  seen  in  Fig.  6. 

5.  Conclusions 

(1)  The  thrust  cannot  be  understated  in  the  transfer  from  a 
P^tking  orbit  to  a  synchronous  transfer  orbit  so  as  to  avoid  an 
overstated  velocity  loss.  If  a  small  thrust  is  to  be  adopted,  a 
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more  complex  control  program  is  required; 

(2)  During  a  transfer  between  non-intersecting  orbits,  a 
small  thrust  propulsion  is  preferable  for  a  primaiy  orbital 
insertion  so  as  to  avoid  too  much  velocity  loss.  In  this  case,  a 
secondary  or  multiple  orbital  transfer  is  generally  adopted. 
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Measurement  of  Spacecraft  Position  and  Attitude 
in  Rendezvous  Docking  Using  Computer  Vision 


Chai  Xiping,  Dai  Yongjiang,  Zhao  Yuan,  Wang  Yan 

(Department  of  Applied  Physics,  Heilongjiang  Engineering 
University,  Harbin,  150002) 


Abstract;  This  paper  introduces  a  method  of  measuring 
spacecraft  position  and  attitude  in  rendezvous  docking. 
With  small  computation  load,  this  simple  method  can 
provide,  through  measurement,  the  position  and  attitude 
angle  of  a  tracking  aircraft  relative  to  a  target 
aircraft. 

Key  words;  computer  vision,  rendezvous  docking, 
relative  position,  attitude  angle 

1 .  Principle  of  Measurement 


To  successfully  accomplish  a  rendezvous  docking,  it  is 
necessary  to  determine  the  position  and  attitude  of  one  aircraft 
relative  to  the  other. 


In  1981,  an  American  named  Fukui  proposed  a  method  of 
measuring  robot  position  and  attitude  using  computer  vision 
technology  in  solving  the  robot  navigation  problem  [2].  He 
adopted  a  standard  rhombic  pattern  as  shown  in  Fig.  1.  Based  on 
this  concept,  a  new  method  and  related  auxiliary  docking 
equipment  used  to  measure  the  relative  position  and  attitude 
angle  between  two  aircraft  are  described  in  this  paper. 


c 


Fig .  1 .  Standard  Pattern  Used  by  Fukui 


The  auxiliary  docking  equipment  shown  in  Fig.  2  is  made  of 
four  semiconductor  lasers  arranged  in  the  fom  of  a  square,  where 
a  three-dimensional  coordinate  system  takes  the  center  of  the 
square  as  its  original  point  and  two  diagonal  lines  respectively 
as  Y  and  Z  axes,  and  a  line  perpendicular  to  this  plane  and 
pointing  to  a  camera  as  X  axis . 


I 

I  X 


/ 

z 

Fig.  2.  Auxiliary  Docking  Equipment 
and  Docking  Coordinate  System 


In  this  coordinate  system,  two  planes,  Pj  and  P2  are 
defined:  one  plane  passes  through  two  points  C  and  D  and  the 
center  of  the  camera,  while  the  other  passes  through  two  points  A 
and  B  and  the  center  of  the  camera  as  shown  in  Fig.  3.  In  this 
way,  a  complex  three-dimensional  problem  can  be  transformed  into 
a  simple  two-dimensional  problem.  Two  visual  angles  are  defined: 
one  is  ff,  a  tension  angle  of  CD  to  CCD  center,  while  the  other  is 


0,  another  tension  angle  of  AB  to  L  as  shown  in  Fig.  4.  In 
addition,  two  polar  angles  are  also  defined:  one  is  p,j3/p,  an 
included  angle  between  the  vertical  line  AB  that  passes  through 
the  original  coordinate  point  and  the  radius  vector  R,  while  the 
other  is  (3,  an  included  angle  between  the  vertical  line  CD  that 
passes  through  the  original  coordinate  angle  and  the  radius 
vector  R  as  shown  in  Fig .  3 .  Thus ,  the  planes  Pi  and  P2 
respectively  contain  a  polar  angle  and  a  visual  angle.  Suppose 
the  projections  of  C  and  D  on  CCD  are  C'  and  D',  then  the 
included  angle  between  LC'  and  LD'  is  also  ♦.  Thus,  If  can  be 
calculated  based  on  simple  geometric  relations.  Similarly,  0 
can  also  be  calculated. 


Fig.  3.  Transformation  of  Three-dimensional 
Problem  to  Two-dimensional  Problem 
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Fig.  4.  Visual  Angles  Defined 
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Fig.  5.  Two-dimensional  Geometric 
Relations  in  Plane  P2 


The  geometric  relations  in  plane  P2  are  shown  in  Fig.  5. 

The  radius  vector  R  divides  AABL  into  two  triangles  Tj  and 
T2,  in  which 


r 


r\  =  2{R-  ^ 


sina  = 


riTjsind^ 

2RW 


The  area  of  Aabl  is  the  sum  of  Tj  and  T2  areas: 


and 


riT.sin^?  =  r,Rsin<7|  4-  r,Rsin^?; 
=  rf  +  r;  -f  2r^r-,cosd 


By  solving  the  above  four  equations, 


cos, a 


(R-  -  W^)xgd 
2RW 


(1) 

(2) 

(3) 

(4) 

(5) 


Similarly,  in  Pi  plane: 


cos/9  = 


2RW 


To  solve  0  and 


(6) 


1*7 


In  Fig.  6,  and  Wj  on  the  image  plane  respectively  are 
projections  of  AM  and  BM  on  the  image  plane.  The  following  can 
be  derived: 


tgd  =  tg(^i  +  di) 


f(W\  4-  W,) 
f  -  IV, Wa 


(7) 


and  can  be  solved  in  the  seime  way. 


A 


Fig.  6.  Relationship  Between  the 
Projection  of  AB  on  Image  Plane  and  0 


Fig .  7 .  Relationship  between  Two  Patterns 


2 .  Equation  Solution 


When  0  and  if  are  solved,  there  are  still  three  unknovm 
quantities  in  Eqs.  (5)  and  (6),  namely  p,  3  and  R.  Even  if  R  is 
known,  P  and  p  still  may  have  double  solutions,  one  positive  and 
one  negative,  which  depends  on  whether  CCD  is  located  to  the  left 
or  right,  above  or  below  the  target. 

To  solve  the  double  solution  problem,  four  more  lasers  of 
the  same  kind  are  added,  which  are  arranged  in  a  square  foinri  just 
the  way  they  are  in  the  auxiliary  docking  equipment  as  described 
above.  The  center  of  the  square  is  positioned  on  Z  axis  at  a 
fixed  distance  from  the  original  coordinate  point  as  shown  in 
Fig.  7.  Subsequently,  the  following  six  equations  can  be 
obtained ; 


,  (R\  - 

cos^i  - 

(8) 

^  {Rl  — 

cosft  -  2R,W 

(9) 

(R--W^)tge, 
cosp,  -  2R,W 

(10) 

(R?  -  vr-)tge. 

COS^Oj  - 

(11) 

Z\Z  =  R.sin/S,  —  Rjsin/?! 
/?iC0S/9,  =  RiCOSji, 


(12) 

(13) 


Here  AZ  is  the  distance  between  two  squares  in  Fig.  7.  To 
solve  this  non-linear  coupled  equation,  six  new  functions 

fj, . .  fg  are  defined  first,  which  respectively  are  the 

differences  between  the  left  sides  and  the  right  sides  of  this 
coupled  equation.  Considering  the  effect  of  noise  in  practical 


“use,  a  non-linear  least  square  method  is  adopted.  The  Jacobi 
matrix  in  the  coupled  equation  is  solved  as 


Fig.  8.  Flow  Chart  of  Computer  Simulation 


Key:  (1)  Initial  input  value 

(2)  Calculate  . fg 

(3)  Output  result 


Fig.  8  is  a  flow  chart  of  computer  simulation,  where  N  is 
the  given  maximvim  iteration  times;  W,  ifz,  0i,  02/  are  all 

measurable;  ^2(1)/  Pi(l)/  P2(l)/  ^2(1) 

respectively  are  given  initial  values  of  quantities  to  be 

measured.  If  the  maximum  values  of  fi, . fg  are  smaller  than 

or  equal  to  10‘®  after  iteration,  then  it  can  be  believed 
approximately  that  the  0i(2),  ^2(2)/  Pi(2)/  P2(2)/  Ri(2)  and 

R2(2)  obtained  now  are  just  Pj,  pi,  p2/  Ri  and  R2  to  be  solved. 


For  instance, 
fi=0.0291 
+2=0.0314 
e,  =0.0321 
^.=0.03229 
W=0.1500 
AZ=0.5000 
Pi(l)=0.0785 
^jl)=0.0785 
^,<1)=0.0785 
jO,  (1)=0.0785 
Ri(l)=10.0000 
R2(1)=10.0000 

After  iteration,  the  final  results  are: 

)3i=0 .663 
j32=0.613 
P,  =0.515 
=0.536 
Ri=8.124 

R2=7 .826 

In  fact,  Pi  and  Pi  are  the  pitch  angle  and  side  slip  angle 
of  the  target  relative  to  CCD,  from  which  the  relative  position 
of  the  target  is  derived  as: 

j  =  R, cos/3, cos/>, , 

y  =  R.cos/5,sin(0, , 

-  =  R,sin^, 


3 .  Conclusions 


This  paper  describes  a  method  of  measuring  the  position  and 
attitude  angle  of  two  spacecraft  by  transforming  a  three- 
dimensional  problem  into  a  two-dimensional  problem.  Based  on 
simple  geometric  relations,  two  coupled  equations  can  be  derived, 
which  contain  three  unknown  numbers .  To  solve  these  coupled 
equations,  one  more  auxiliary  docking  equipment  of  the  same  kind 
is  added.  The  required  information  about  the  position  and 
attitude  angle  can  be  acquired  after  iteration. 

This  paper,  a  project  funded  by  State  863  Program,  was  received 
on  May  27,  1993. 
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